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In this work we present the resuhs of theoretical analysis of the de Haas- van Alphen oscillations in 
quasi-two-dimensional conductors. We have been studying the effect of the Fermi-liquid correlations 
of charge carriers on the above oscillations. It was shown that at reasonably low temperatures and 
weak electron scattering the Fermi-liquid interactions may cause noticeable changes in both ampli- 
tude and shape of the oscillations even at realistically small values of the Fermi-liquid parameters. 
Also, we show that the Fermi-liquid interactions in the system of the charge carriers may cause 
magnetic instability of a quasi-two-dimensional conductor near the peaks of quantum oscillations 
in the electron density of states at the Fermi surface, indicating the possibility for the diamagnetic 
phase transition within the relevant ranges of the applied magnetic fields. 
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I. INTRODUCTION 

Magnetic quantum oscillations P, 0, Q have been rec- 
ognized as one of the major tools to map Fermi surfaces 
(FS) in metals. Analysis of the experimental data is 
now a well established procedure and is based on the 
classical paper by Lifshitz and Kosevich (LK) ^ who 
first layed out a quantitative theory of the de Haas-van 
Alphen effect. The LK expression for the oscillating part 
of the thermodynamic potential and magnetization of 
metallic electrons in a strong (quantizing) magnetic field 
B was derived assuming that conduction electrons are 
noninteracting quasiparticles in a periodic crystal poten- 
tial. This potential determines the electron dispersion, 
_B(p) (p being the electron quasimomentum) , and there- 
fore, the effective mass of conduction electrons m* and 
their FS. In fact, conduction electrons interact with each 
other. Studies of modifications of the LK results arising 
due to electron-electron interactions within the general 
many-body quantum field theoretical approach started 
in early sixties in the works of Luttinger [HI, and con- 
tinued through the next three decades [o]. It was shown 
that electron-electron interactions may bring noticeable 
changes in the de Haas-van Alphen oscillations which 
makes further analysis worthwhile. 

One of the oldest and still powerful method to deal 
with electron-electron interaction is the Landau Fermi- 
liquid (FL) theory It is important to realize that 
while the phenomenological Fermi liquid theory and the 
microscopic many-body perturbation theory (and, when 
applicable, the exact density functional theory) by def- 
inition lead to the same observable quantities, such as 
response to an external field, both zero approximation 
and its renormalization depend on the taken approach. 
A discussion to this effect in application to the dielectric 
response of metals can be found, for instance, in Ref. [13] ■ 
It is always instructive to look at the same phenomenon 



from different points of view. Respecting the value of 
the many-body perturbation theory as applied to quan- 
tum oscillations 6], we emphasize that the Fermi- liquid 
theory has provided important insights in such areas, rel- 
evant for the de Haas- van Alphen physics, as high fre- 
quency collective modes in metals [HI, [l2, 13, 14, 15, 1^ . 
or oscillations of various thermod yna mic observables in 
quantizing magnetic fields 17, 18, 3, 2ol |. 

An advantage of this phenomenological theory is that 
it enables to describe the effects of quasiparticles in- 
teractions in such a way that makes the interpretation 
of the results rather transparent, as compared with the 
field-theoretical methods. At the same time the many- 
body approach brings general but cumbersome results, 
and usually it takes great calculational efforts and/or 
significant simplifications to get suitable expressions for 
comparison with experimental data. Which is more im- 
portant, adopted simplifications may lead to omission of 
some qualitative effects of electron-electron interactions 
in quantum oscillations, as we show below. 

In the last two decades an entire series of quasi-two- 
dimensional (Q2D) materials with metallic type conduc- 
tivity has been synthesized. These are organic conduc- 
tors belonging to the family of tetrathiafulvalene salts, 
dichalcogenides of transition metals, intercalated com- 
pounds and some other. At present, these materials at- 
tract a significant interest. Their electronic properties are 
intensively studied, and the de Haas- van Alphen effect is 
employed as a tool in these studies [1, [1] . Correspond- 
ingly, the theory of this effect in the Q2D materials is 
currently being developed [HI, HI, IHlIi, IHIiil . The 
present analysis of the effect of Fermi-liquid interactions 
on the de Haas-van Alphen oscillations contributes to 
the above theory. Also, the analysis is motivated by the 
special features in the electron spectra in Q2D materials 
providing better opportunities for the Fermi-liquid effects 
to be manifested, as was mentioned in some earlier works 
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In this paper we show how renormalizations of con- 
duction electron characteristics arising from FL interac- 
tions affect quantum oscillations, and we express them in 
the form appropriate for comparison with experiments. 
Also, we show that a magnetic phase transition lead- 
ing to emergence of diamagnetic domains may happen 
at low temperatures when the cyclotron quantum fku 
is very large compared to the temperature expressed in 
the energy units: {huj ^ ksT) . We emphasize that 
the analyzed effects are different from the many-body 
renormalizations of the band structure. Within the phe- 
nomcnological theory the latter are already included in 
the ground state of conduction electrons. 



II. FERMI-LIQUID RENORMALIZATIONS OF 
THE CONDUCTION ELECTRON 
CHARACTERISTICS 

Within the phenomenological Landau FL theory single 
quasiparticle energies get renormalized, and the renor- 
malization is determined by the distribution of excited 
quasiparticles. Accordingly, the energy of a "bare" (non- 
interacting) quasiparticle £'o(p) moving in the effective 
crystal potential is replaced by the renormalized energy 
defined by the relation Q: 

E^ip, r,t) = £;o(p) + J2 ^(P' P'' ^')^PiP'^ s'. r.O- (1) 

p'cr' 

Note that Eq{p) here is the energy spectrum in the ab- 
sence of any excited quasiparticles, that is, at equilib- 
rium and at zero temperature, while Sp{p' ,s' ,r,t) rep- 
resents the nonequilibrium part of the electrons distri- 
bution function, which may depend on both position r 
of the quasiparticle and time t. Also, s , s' are spin 
Pauli matrices, (cr is the spin quantum number), and 
F(p,s;p',s') is the Fermi-liquid kernel (Landau corre- 
lation function) , which describes additional renormaliza- 
tion of the quasiparticle spectrum due to interaction with 
other excited quasiparticles (but not with all electrons 
in the system, which is included in Eo{p) ). Neglecting 
spin-orbit interactions, the Landau correlation function 
may be written as: 

F(p, s; p', s') = ^(p, p') + 4ij{p, p')(ss'). (2) 

As follows from the Eq. [H the conduction electron 
velocity v — WpE differs from the bare velocity vq = 
Vpi^Q. To proceed in our analysis we need to bring ve- 
locities V and Vq into correlation. For brevity we do 
not explicitly write out the variables r,t in the Eq. [1] 
in further calculations. This omission does not influence 
the results. Differentiating the Eq. [l]we obtain: 



v^(p) = Vo(p) + Vp 



■F{p,s;p',s')Sp{p\s') 



pV 



(3) 



The electron distribution function p(p, s) is the sum of 
the equilibrium part pQ (the latter coincides with the 
Fermi distrubution function for quasiparticles with single 
particle energies Eo{p)) and the nonequlibrium correc- 
tion Sp. Multiplying both parts of the Eq. [3] by p(p, s) 
and performing summation over p, a we get: 

^p(p, s)v„(p) = p{p, s)vo(p) + pip, s)Vp 



per 



per 



per 



^F(p,s;p',s')Mp',s') 

p'cr' 



(4) 



The second term on the right hand side of the Eq. [4] 
could be converted to the form 



E 

per 



^F(p,s;p',s')Mp',s') 

p'cr' 



Vpp(p,s). (5) 



Keeping only the terms linear in Sp (which is sup- 
posed to be small compared to the equilibrium part of 
the distribution function), we may approximate Vp/5 as 
Vo-(p) f^'"^ . Here, / is the Fermi distribution function 



'9-Bp 

and the quasiparticle energies Ea-{p) correspond to the 
local equilibrium of the electron liquid. Also, assuming 
that the FS of a considered metal possesses a center of 
symmetry, we get 

pa pa 

E P^P' s)vo(p) = E '^^(P' ^)vo(p)- (6) 

per per 

Then, using the well known relation F(p, s;p',s') = 
F(p',s';p,s) and carrying out the replacement 
p, s ^ p', s' in the sums included in the Eq. [51 we could 
rewrite Eq. 2] as 



per 



^-(p) + E l^^(P' p'' (p') 

per ' 



(7) 



Solving this for v , we obtain 
v.(p)-Vo(p)- ^ ||^i^F(p,s;p',s')v.,(p')- (8) 

p',er' P 

Exact expressions for the functions (p{p, p') and 
'0(p, p') are of course unknown. The simplest approx- 
imation is to treat them as constants. This approxima- 
tion is reasonable as long as the interaction of quasipar- 
ticles (located at r and r', respectively), is extremely 
short range, so that the interaction can be approximated 
as V{r,r') = IS{r — r'). Using this approximation one 
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captures some FL effects but in general case it is not 
sufficient. 

As a next step, one may expand tlie Fermi-liquid func- 
tions in the Eq. [2]in basis functions respecting the crystal 
symmetry, such as Allen's Fermi surface harmonics [27| : 

d dj 

j — l m— 1 

d dj 

i^ip, p') = E E i^,{p.p')R,m{0, mUO', *')• (9) 

j — l m—1 

Here, we introduce spherical coordinates for p : p = 
{p, 0, $); d is the order of the point group; index j labels 
irreducible representations of the group; dj is the dimen- 
sion of the j -th irreducible representation; {Rjm{d,^)} 
is a basis of the j -th irreducible representation including 
dj functions. 

For an isotropic metal the spherical harmonics Yjm 
can be used as the basis. Including orbital moments up to 
j = 2 we have, for a cubic symmetry (cubic harmonics): 



{PzP'zPxP'x +PzP'zPyP'y +PxP'xPyP'y) 



HpI-pIW^-p'') 



+ i'^pl ~vl~ pI){2p'I - P'^ - P'^). (10) 

The coefficients ip^ip are material dependent constants. 
A common feature of Q2D metals is their layered struc- 
ture with a pronounced anisotropy of the electrical con- 
ductivity. In such materials electron energy only weakly 
depends on the quasimomentum projection p = pn on 
the normal n to the layers plane. In further considera- 
tion we assume n = (0, 0, 1) and we neglect the asymme- 
tries of the electron spectrum in the layers planes. Then 
the relevant Fermi surface is axially symmetrical. 

For systems with an axial symmetry this expression 
(fTU)) needs to be correspondingly modified. For instance, 
in the first order we have 



V(P,P') 
V'(P,P') 



<y5io 



PzPz 



(fill 
Ipll 



{PxP'x+PyP'y)- 



(11) 



This expression will be used from now on in the present 
paper. 

When an external magnetic field B = (0, 0, B) is ap- 
plied the spin degeneracy of the single electron energies 
is lifted, and we can write: 



where Eq (p) does not depend on the electron spin, g is 
the electron Lande factor, and /3o — eh/2mQC is the Bohr 
magneton (mo is the free electron mass). The nonequi- 
librium correction to the electron distribution function 
satisfies the equation [t^: 



(5p(p,s) = (5p(p,s) 



5/, 



pa 



^i^(p,s,p',s')Mp',s') 

p',cr' 

_ (13) 
where i5/9(p,s) describes the deviation of the electron 
liquid from the state of local equilibrium. When the de- 
viation arises due to the effect of the applied magnetic 
field = -{dfpa/dEp„)g(iaB. Substituting Eq. [U] 
into Eq. [13] and using the result in Eq. [l]we get: 



AE = AEq - b*agl3oB 



l + 6n ■ 



(14) 



where h* = 6o/(l -I- 6o)i ^nd feg is a dimensionless 
parameter describing FL interactions of the conduction 
electrons, namely: fog = — i'o(0)^/'o where t'o(O) is the 
density of states of noninteracting conduction electrons 
on the Fermi surface in the absence of the magnetic 
field. This is nothing but the standard Stoner renor- 
malization of the paramagnetic susceptibility. Note that 
here tpQ plays the role of the Stoner parameter / — 
(S'^Exc/Sp^Spi) in the density functional theory, or of 
the contact Coulomb interaction in the many-body the- 
ory. 

The Luttinger theorem dictates the Fermi surface vol- 
ume. Therefore, the radius and the cross-sectional areas 
of the Fermi sphere associated with an isotropic Fermi 
liquid remain unchanged due to quasiparticles interac- 
tions. In realistic metals whose conduction electrons 
form anisotropic Fermi liquids one may expect some mi- 
nor changes in the FS geometry to appear. Such effects 
could be considered elsewhere. In the present work we ne- 
glect them. So, in further consideration we assume that 
FL interactions do not affect the FS geometry. Then 
the cross sectional areas of the Fermi surface A{pz) cut 
out by the planes perpendicular to the magnetic field B 
do not change when electron-electron interactions are ac- 
counted for. However, the cyclotron masses of conduction 
electrons undergo renormalization due to the electron- 
electron interactions. The cyclotron mass is defined as: 



1 OA 
'^^=2^dE 



dl 



(15) 



Here, dl = ^ dp^ + dpy is the element of length along 
the cyclotron orbit in the quasimomentum space, v± = 



EoAp) - Eo{p) + agPoB = Eo{p) + AEq 



(12) 



+ Vy, Va = dE/dpa (ck = x,y), and /i is the chem- 
ical potential of conduction electrons. 

Substituting Eqs. [TT] into the Eq. [5] we get v± = 

w^o/(l+ai) where Vj_o = ^v^q + u^q, and ai is related 
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to the FL parameter ipu as follows: 
- i^o{0)pl(pii/3 ^ ai 



(16) 



where po is the maximum value of the longitudinal com- 
ponent of quasimomcntum. So we get: 



mj_ = ?7i±o(l + ai), 



(17) 



m±Q being the cyclotron mass of noninteracting quasi- 
particles. In the case of isotropic electron system the 
cyclotron mass to^o coincides with the crystalline ef- 
fective mass m* Therefore, our result agrees with the 
standard isotropic FL theory. 

Other quantities, such as the chemical potential of con- 
duction electrons and their compressibility, may expe- 
rience different renormalizations, as well. The latter, 
for instance, is renormalized by a factor 1/(1 + ao) = 
1/[1 — vo{0)^o] , and the former by the factor (1 -I- ao) 
[7|. So, the renormalized density of states z^(0) appears 
in the expressions for the electron compressibility and the 
velocity of sound in metals. 

The model of the extremely short range (contact) 
Coulomb interaction between quasiparticles is often em- 
ployed, while applying the many-body theoretical ap- 
proach to study de Haas-van Alphen effect (see e.g. Ref. 
[a|)- Within the phenomenological FL theory this model 
results in the approximation of the functions (p(jp, p') 
and ?/'(p,p') by constants ipo and ipo, respectively. 
Such approximation enables us to get the Stoner renor- 
malization of the paramagnetic susceptibility, and elec- 
tron compressibility as shown above. However, it misses 
Fermi-liquid effects associated with the subsequent FL 
coefficients included in the Eqs. [SHTTl which could be 
significant in renormalizations of other parameters char- 
acterizing the charge carriers such as their cyclotron 
masses. 



III. QUANTUM OSCILLATIONS OF THE 
LONGITUDINAL VELOCITY OF 
CHARGE-CARRIERS IN Q2D CONDUCTORS 

In further calculations wc adopt the commonly used 
tight-binding approximation for the charge carriers spec- 
trum in a quasi-two-dimensional metal. So, when a quan- 
tizing magnetic field is applied, the charge carriers ener- 
gies may be written in the form: 



-\+(jhwn-2tcofi (tt— ] . (18) 
2/ V Po, 



where Hluq is the spin splitting energy, t is the inter- 
layer transport integral, and po — nh/ L where L is the 
interlayer distance. This expression (18) describes single 
particle energies of noninteracting quasiparticles. Now, 
the relation of matrix elements of renormalized v^^' and 



bare Voi/,y' velocities in accordance with Eq. [8] takes on 
the form 91: 



(19) 



Here, E^ is the quasiparticle energy including the correc- 
tion arising due to the FL interactions, and v = {a, a} 
is the set of quantum numbers of an electron in the mag- 
netic field. The subset a includes the orbital num- 
bers n^Pz and a;o (the latter labels the positions of 
the cyclotron orbits centers). Also, F^^,"^ — if^l^^ + 
4-!/'"^?^ (sSj^) are the matrix elements of the Fermi- liquid 
kernel. 

For an axially symmetrical FS the off-diagonal matrix 
elements of the longitudinal velocity vanish and we ob- 
tain: 



(20) 



Substituting the expression for the Fermi-liquid kernel 
into the Eq. [201 we get: 



(21) 



where both Vaa and v^^ only depend on pz, so in 
the further calculations we will use the notation Vaa = 
v{pz), = v''{pz)- These matrix elements could be 
found from the system of equations that results from Eqs. 
[201 and mj 



(22) 



(23) 



Here, 
r 



<Ji. (24) 



The de Haas- van Alphen oscillations are observed in mag- 
netic fields when the Landau levels spacing is small com- 
pared to the chemical potential of electrons {hui <C /i). 
Under these conditions we may approximate the Fermi- 
liquid kernel by its expression in the absence of the mag- 
netic fields (Eq. [2]). Using the above-described approx- 
imation of the Fermi-liquid functions (p{p,p') ■0(p,p') , 
Eqs. [Til we can solve Eqs. [221 [23 To this end, we need 
some averages over the Fermi surface, namely: 



1 



R^-J2^^^MPz)Pz - ^^^^2 

a 

df{E,,^APz)) 



E 



dEn,a{Pz 



-VQiPz)Pzdpz 



(25) 
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(26) 



The Fermi- liquid effects enter the system [221I231 through 
the averages A,A',B,B' closely related to the Fermi- 
liquid parameters: 



otaPz 1 



^'^-^loEToaP'- (27) 



The expressions for B, B' could be obtained replacing 
<Pio by Vio in Eqs. W\ 

Applying the Poisson summation formula, 

(p{n) — / exp(27rir?i)iy9(n)(iri. (28) 



to Eq. [H we get: 



R = - 



CO ^ 

1 + 2Re^ exp(27rirn) I . 



Vo{pz)Pz 



(29) 



So, we see that oscillating terms appear in the expres- 
sions for R and other averages over the Fermi surface 
included in Egs. [22l[23l Due to this reason, an oscillating 
term occurs in the resulting formula for the renormalized 
longitudinal velocity Va-{pz)- This oscillating term origi- 
nates from the Fermi-liquid interactions between charge 
carriers, and it appears only when the FL coefficients 
ifio and "010 are taken into account, that is, beyond the 
contact approximation for the Coulomb interaction. 

We get the following results for the oscillating parts of 
i?, R' : 



i? = 2iV ( I 1 A, 



R'^2N{j] A^ (30) 



where N is the electron density, and the functions A 
and A'* have the form: 



A = V ^—^D{r) sin ( 27rr^ ) cos ( nr— 



r=l 

X Jo ( Anr 



B 



huj* J ' 



A« =Y_^L2]LD{r) cos (^Svrr^^ sin (^Trr^ 



(31) 



X Jo Anr 



t 

huj* 



(32) 



Here, F ~ cAo/2Trhe, Aq is the FS cross-sectional area 
at Pz = ±po/2; and the cyclotron quantum ftw* and 
spin-splitting energy Hloq are renormalized according to 
Eqs. [HI [T71 The damping factor D{r) describes the 
effects of the temperature and electron scattering on the 



magnetic quantum oscillations, and Jo (a;) is the Bessel 
function. The simplest and well known approximation 
for D{r) equates it to the product RT{r)Rr{r) where 
Rxir) — rxj sinh(rx) {x — 2Ti'^kBT /hw*) is the temper- 
ature factor and i?T(r) — exp[— 7rr/w*r] is the Dingle 
factor describing the effects of electrons scattering char- 
acterized by the scattering time r. The temperature fac- 
tor appears in the Eqs. I31I32I as a result of standard 
calculations repeatedly described in the relevant works 
starting from the LK paper [4]. The Dingle factor can- 
not be straightforwardly computed starting from the ex- 
pessions like [25j [26l This term is phenomenologically 
included in the Eqs I31|32l in the same way as in the 
Shenberg's book Under low temperatures required 
to observe magnetic quantum oscillations, the value of r 
is mostly determined by the impurity scattering. 

Using the microscopic many-body perturbation the- 
ory it was shown that in this case the Dingle term re- 
tains its form, and the corresponding relaxation time 
could be expressed in terms of the electron self-energy 
part S arising due to the presence of impurities, namely: 
T^^ = 21mTi/h. In strong magnetic fields the self-energy 
E gains an oscillating term which describes quantum os- 
cillations of this quantity [H, H^. So, the scattering 
time becomes dependent of the magnetic field B. A thor- 
ough analysis carried out in the earlier works of Cham- 
pel and Mineev [2^ and Grigoriev [1^ shows that the 
oscillating correction to the scattering time could be ne- 
glected when the FS of a Q2D metal is noticeably warped 
(47rt > hw*). In such cases one may treat t as a phe- 
nomenological constant. However, when the FS is very 
close to a pure cylinder (47rt ^ hio*) the scattering time 
oscillations must be taken into consideration in studies 
of the de Haas- van Alphen effect. These oscillations may 
bring some changes in both shape and magnitude of the 
magnetization oscillations but we do not discuss the issue 
in the present work. In further analysis we assume that 
Airt > Huj*. 

One may notice that the oscillating function A has 
exactly the same form as that describing the magneti- 
zation oscillations in Q2D metals when the Fermi-liquid 
effects are omitted from the consideration (see e.g. Refs. 
23l . 24 1). Also, the Fermi-liquid terms included in the ex- 
pression (P?]) exhibit oscillations in the strong magnetic 
field. For instance, applying the Poisson summation for- 
mula to the expressions (j27|) we can convert these expres- 
sions to the form: A = ai{l + S), A' = aiS'' where the 
oscillating functions 6 and (5* are 



..|(-:,.-z„.-)».(2„|)c„.(„^)5(i=l) 



(33) 



6 



(34) 



The factors 5 and Q entered in Eqs. [321 and | 
expressed in the series of the Bessel functions: 



are 



m— 1 



^-^J2M, (35) 



^2 (2m +1)2 

m— 



J; 



2m+l 



(x). 



(36) 



The expressions for B, B' arc similar to those for A, A' 
and we may get to former by replacing the factor ai by 
another constant bi . The oscillating function S behaves 
like the function describing quantum oscillations of the 
charge carriers density of states (DOS) on the FS of a 
Q2D metal (see Appendix). As for the parameters ai,bi 
we can define oi = — fo(0)poi^io/3 and bi is similarly 
defined, namely: bi = — i'o(0)poV'io/3- We remark that 
the parameter oi differs from ai which enters the ex- 
pression for the cyclotron mass (see Eq. [17)) . This reflects 
the anisotropy of electron properties in Q2D conductors. 



IV. QUANTUM OSCILLATIONS IN THE 
MAGNETIZATION 



To compute the longitudinal magnetization Af| 
start from the standard expression: 



M\\{B,T,^i) = ^u{B,T,^i)^-[—^ 



we 



(37) 



Here, the magnetization depends on the temperature T 
and on the chemical potential of the charge carriers /z, 
and H is the external magnetic field related to the field 
B inside the metal as B = H + 47rM. When the mag- 
netic field is directed along a symmetry axis of a high 
order we may assume that the fields B and H are par- 
allel. One may neglect the difference between B and 
H when the magnetization is weak. Otherwise the uni- 
form magnetic state becomes unstable, and the Condon 
diamagnetic domains form, with the alternating signs of 
the longitudinal magnetization [i^. We will discuss this 
possibility later. Now, we assume by B in the Eq. 
[37l To incorporate the effects of electron interactions we 
assume, in the spirit of the PL theory, that the thermo- 
dynamic potential f2 has the same form as for noninter- 
acting quasiparticles, but with the quasiparticle energies 
fully rcnormalizcd by their interaction: 



-/cfiT^lni 1 + cxp 



knT 



(38) 



In this expression E^, is the quasiparticle energy includ- 
ing the correction arising due to the PL interactions, and 
fcs is the Boltzmann's constant. 

Accordingly, we rewrite Eq. [38las follows: 



/ In < 1 + cxp 



dpz 
(39) 

where — hc/eB is the squared magnetic length. Per- 
forming integration by parts, Eq. 1391 becomes: 



n = - 



(40) 



Applying the Poisson summation formula, we get: 



^3^Y1 J J '^P^f(^rL,aiPz))VaiPz)Pz 



l + 2Re^exp[2 7 



(41) 



r=l 



So we see that the expression for the thermodynamic 
potential includes two oscillating terms. One originates 
from the oscillating part of v„{pz). The second term in- 
side the braces in the Eq. [41] gives another oscillating 
contribution. 

The effects of temperature and spin splitting on the 
magnetic oscillations are already accounted for in the Eq. 

Assuming Ant > huj, we take into account the effect 
of electron scattering adding an imaginary part ih/2T to 
the electron energies fll|. After standard manipulations, 
we obtain the following expression for the oscillating part 
of the longitudinal magnetization: 



AM|| =-2A^/3^(l-3aD 
II ^* 

A - 3{al + bl)S + 3bl{AS - A'S") - 9albl{S^ ~ S"^) 



1 + 3{al + bl)S + 9a*5*((52 - ^-2) 



(42) 



where = ai/(l + 3ai); bl = 6i/(l + 36i). This is 
the main result of the present work. It shows that the 
Fermi-liquid interactions may bring significant changes 
in the de Haas-van Alphen oscillations. Below, we an- 
alyze these changes. If we may neglect the oscillating 
corrections proportional to , 6* , then our result for 
AAf|| reduces to the usual LK form with some renor- 
malizations arising from the quasiparticle interactions. 
The cyclotron mass m± differs from the bare Fermi liq- 
uid cyclotron mass before the quasiparticle interaction is 
taken into account (cf. Eq. [T7)) . and Hujq includes the 
extra factor (l-t-6o)~^- Also, the factor {l — 3al) mod- 
ifies the magnetic oscillations magnitudes. As for the 
oscillations frequencies, they remain unchanged by the 
PL interactions, as expected. 
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V. DISCUSSION 



Comparing our result ([32]) with the corresponding re- 
sult reported by Wasserman and Springfield 6], we see 
that these results agree with each other. A seeming dif- 
ference in the expressions for the oscillations frequencies 
arises due to the fact that in Ref. the frequencies 
are expressed in terms of the chemical potential of elec- 
trons fi instead of the cross-sectional areas of the Fermi 
surface. It is worth reiterating that "unrenormalized" 
mass in the FL theory is already renormalized (sometimes 
strongly) from fully noninteracting (or density functional 
- calculated) mass. Again, we remark that the present 
analysis was carried out assuming noticeable/significand 
FS warping {Ant > huj*), so, we may neglect the mag- 
netic field dependence of the electrons scattering r treat- 
ing the latter a constant phenomenological parameter. 
This results in a simple form of the Dingle damping fac- 
tor Rrir) describing the effects of electrons impurity 
scattering. 

The LK form of the expression for the longitudi- 
nal magnetization is suitable to describe de Haas-van 
Alphen oscillations in conventional three dimensional 
metals within the whole range of temperatures. However, 
this is not true for quasi-two-dimensional conductors. 
The Fermi surface of such a conductor is nearly cylindri- 
cal in shape, therefore the oscillating term in the denom- 
inator of Eq. |42] significantly increases. The oscillations 
of the denominator of Eq. [42] occur due to the functions 
6 and 5'^. These functions are presented in the Fig. 1, 
and we see that at low temperatures and weak scattering 
\n{hoj*/kBT*) > t/hw* (T* + Td, Td = h/2TrkBT 
is the Dingle temperature) the peak values may be of 
the order of 1 , especially for a rather weakly warped 
FS {t/huj* ^ 0.1 0.5). So, quantum oscillations in the 
magnetization in the electron Fermi-liquid in quasi-two- 
dimensional metals may have more complicated struc- 
ture than those in the electron gas described by the 
LK formula The effect of the Fermi-liquid interac- 
tions on these oscillations depends on the values of the 
Fermi-liquid parameters a* , and on the damping fac- 
tor D{r) = RT{r)Rr{r) included in the expressions for 
the oscillating functions. The FS shape determined by 
the ratio t/huj is important as well. 

The most favorable conditions for the changes in the 
magnetization oscillations to be revealed occur when the 
oscillating terms in the denominator of the Eq. (42) may 
take on values of the order of unity at the peaks of oscil- 
lations. We may estimate the peak values 5m and (5^^ of 
the functions included into the above denominator using 
the Euler-Macloren formula. The estimations depend on 
the shape of the FS of the Q2D conductor. When the FS 
is significantly crimped (t S> tiuj) we obtain: drmS^ ^ 
{nuj/tf/'^{kBT*)-^/'^. So, we may expect the Fermi- 
liquid interaction to be distinctly manifested in the mag- 




(B(j/B-1)x10 



(Bjj/B-1)x10'^ 



FIG. 1: The magnetic field dependencies of the functions 5 
(solid lines) and 5" (dashed lines). The curves are plotted 
at Bo = lOr, F/Bo = 300, 2T?e* jhw* = 0.5 for f/ficj* = 2 
(left panel) and t/htu* = 0.3 (right panel). 



netization oscillations when |aJ|(Sa;/i)^/^(fcBT*) ~ 
1 or \bl\(huj/tf/'^(kBT*)-^/'^ - 1 or both. In aU prob- 
ability the Fermi-liquid parameters are small in magni- 
tude (|a*|, \h\\ <C 1). Nevertheless, the changes in the de 
Haas-van Alphen oscillations arising due to the Fermi- 
liquid effects may occur at UbT* <^ 1. It is worthwhile 
to remark that due to the character of the electron spec- 
tra, the Q2D conductors provide better opportunities for 
observations of the Fermi-liquid effects in the de Haas- 
van Alphen oscillations than conventional 3D metals. In 
the latter the peak values of the oscillating functions 5, 6^ 
have the order [hw / jiY/'^iJtBT*)^^/'^ . Typical values of 
the transfer integral t are much smaller than those of 
the chemical potential /i, therefore significantly smaller 
values of fcsT* and/or greater values of the parameters 
a{,b\ are required for the Fermi-liquid effects to be re- 
vealed in 3D metals. 

Due to the special character of the electron spectra in 
the Q2D metals, the variations in the magnetization os- 
cillations may be noticeable at reasonably small values of 
the Fermi-liquid constants. In the Fig. 2 we compare the 
oscillations arising in a gas of the charge carriers (top left 
panel) with those influenced by the Fermi-liquid interac- 
tions between them. All curves included in this figure are 
plotted within the limit t > huj. We see that both magni- 
tude and shape of the oscillations noticeably vary due to 
the Fermi-liquid effects. When t/hio* ~ 0.1-^-0.5 the FS 
shape is closer to a perfect unwarped cylinder, the mag- 
netization oscillations accept the well known sawtoothed 
shape, shown in the Fig. 3. Again, when the Fermi-liquid 
interaction produced terms are included into the expres- 
sion for AM, this bring some changes in the magnitude 
and shape of the oscillations. These changes are more 
significant when (a^ + ^i) < 0. 

The most important manifestation of the Fermi-liquid 
effects occurs in very clean conductors at low tempera- 
tures when T* is reduced so much that (a* + bl)Sm is 
greater than 1. Then the denominator of the Eq. [55| 
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FIG. 2: The effect of the Fermi-hquid interactions on the de 
Haas-van Alphen oscillations in a Q2D metal at t/huj* — 2. 
The curves are plotted using Eq. |42l Mo = 2Np. Calcu- 
lations are carried out for al — bl = (top left panel) 
al = bl = 0.02 (top right panel), al = bl = -0.02; -0.04 
(bottom left and right panel, respectively). The remaining 
parameters are the same as used in the figure 1. The dashed 
lines in the top right panel and in the bottom panels represent 
oscillations in the system of noninteracting quasiparticles. 



FIG. 3: The effect of the Fermi- liquid interactions on the de 
Haas-van Alphen oscillations in a Q2D metal at t/hcj* = 0.3. 
The curves are plotted using Eq. |42]at al = bl = (top left 
panel), al — bl — 0.02 (top right panel), al — bl — —0.02 
(bottom left panel), and a* = bi — —0.04 (bottom right 
panel). The remaining parameters are the same as in the 
figure 1. Dashed lines represent the oscillations in the gas of 
charge carriers. 



becomes zero at some points near the peaks of the DOS 
oscillations provided that (a^ + 6*) < 0. This is illus- 
trated in the Fig. 4. Correspondingly, AM diverges at 
these points which indicates the magnetic instability of 
the system. 

This means that the condition for the uniform mag- 
netization of the electron liquid is violated near the os- 
cillations maxima, and the diamagnetic domains could 
emerge. It is known that both crystal anisotropy and de- 
magnetization effects originating from the shape of the 
metal sample, could modify the relation between B and 
H and cause magnetic instability which results in the 
occurence of the diamagnetic domains [28]. Our result 
demonstrates that the interactions of conducting elec- 
trons also may play an important role in magnetic phase 
transitions. In principle, such magnetic instabilities may 
appear in 3D metals as well which was shown earlier ana- 
lyzing quantum oscillations in the longitudinal magnetic 
suscegtibility of the isotropic electron liquid (see Refs. 
[H Q). However, we may hardly expect these tran- 
sitions to appear in conventional metals for the require- 
ments on the temperature and intensity of scattering pro- 
cesses are very strict. Estimations made in the earlier 
works 2^, 2^ show that the temperature (including the 



Dingle correction) must be about lOmK or less for these 
diamagnetic phase transitions to emerge in conventional 
metals. On the contrary, the special (nearly cylindri- 
cal) shape of the FS in Q2D conductors gives grounds to 
expect the above transitions to appear in realistic exper- 
iments. 

To summarize, in the present work we theoretically 
analyzed possible manifestations of the FL interactions 
(that is, residual interactions of excited quasiparticles) in 
the de Haas- van Alphen oscillations in Q2D conductors. 
The same approach can be easily applied for a metal 
with the crystalline lattice of arbitrary symmtery, using 
the appropriate basis for expanding the FL functions. 
So, the phenomenological Fermi-liquid theory becomes 
more realistic and suitable to analyze effects of electron 
interactions in actual metals. 

We showed that the residual quasiparticle interactions 
affect all damping factors inserted in the LK formula 
through the renormalization of the cyclotron mass. The 
spin splitting is renormalized as well, in a manner similar 
to the so-called Stoner enhancement. The frequency of 
the oscillations remains unchanged for it is determined 
with the main geometrical characteristics of the Fermi 
surface, which probably are not affected by electron- 
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FIG. 4; The plot of the function Y = (a^ + &t )(S + atbj: ("^^ - 
5"^) near the peak of the DOS quantum oscillations at 
t/huj* = 0.3 which illustrates that the denominator in the 
Eq. |42]may become zero in the vicinities of these peaks (left 
panel) . The divergencies in the oscillating part of magnetiza- 
tion described by Eq. |42] indicating the magnetic instability 
caused by the Fermi-liquid effects (right panel). The curves 
are plotted assuming a* = b* = —0.04. The remaining pa- 
rameters have the same values as used in the figure 1. 



electron interactions. However, the shape and magni- 
tude of the oscillations are affected due to the FL effects, 
and their changes may be noticeable. Also, the obtained 
results indicate that (under the relevant conditions) the 
electron interactions may break down the magnetic sta- 
bility of the material creating an opportunity for the dia- 
magnetic phase transition. The discussed effects may be 
available for observations in realistic experiments bring- 
ing extra informations concerning electronic properties of 
quasi-two-dimensional metals. 

Finally, we want to emphasize once again that the 
renormalizations, most importantly, mass renormaliza- 
tion, are in addition to what is conventionally called 
"mass renormalizatin" , namely, renormalization of the 
specific heat coefficient compared to band structure cal- 
culations. It is usually implicitely assumed that the 
weighted average of the de Haas-van Alphen mass reno- 
ramization is exactly equal to the specific heat renor- 
malization, i.e., the FL effects are small. In many cases 
this is a good approximation, but one can never exclude 
a possibility that in some materials these two masses 
may be different, namely, the de Haas- van Alphen mass 
may be larger. A curious example when one would have 
needed to exercise caution, but did not, is given by Ref. 
[30| . where quantum oscillations in a highly uncoven- 



tional metal, Naa;Co0 2, were measured, and it was 
taken for granted that the large mass renormalization 
found in the experiment should be fully accounted for 
in specific heat. Based on this assumption, a natural 
and straightforward interpretation of the data was aban- 
doned and counterintuitive explanation, requiring some 
unverified assumptions was accepted. It is possible that 
the results reported in the Ref. [s^ give a case where 
additional mass renormalization discussed in this paper 
is significant. Hopefully, at some point we will see a care- 
ful and accurate experimental study on various materials 
that would compare the de Haas- van Alphen masses with 
the thermodynamic masses and give us a quantitative an- 
swer on how different may the two be in real life. 



APPENDIX 

Here, we analyze the expressions for the oscillating 
functions S, (5* within the limits of significant (t ^ huj) 
FS warping. We may use the standard asymptotics for 
the Bessel functions , at x 3> 1, namely: 



Jkix) 



Trk TT 
T ~ 4 



Substituting these aproximation into Eqs. 
tain: 



S{x) 
Q{x) 
where 
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(43) 
we ob- 

(44) 
(45) 
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m— 1 

So, we have: 
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Q{x) 



4 V TTZ 

3 FY 



cos 



4 V TTX 



sm 



^- (46) 



(47) 
(48) 



J 



Using these results we may write the following expressions for 5,5^ at t ^ hu 



5 f hio* 
~4 [27rH 
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(-ir 



D{r) cos 



4 \2^t 
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1/2 OO 



" F' 
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7rr— ^ 
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B_ 




hw* 
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LU* 



27rr- 



F 
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Anrt n 
hjj* 4 



(49) 
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Or: 



S = 



\ 1/ Z i-X-J 



2itHJ ^ Jr 
Likewise, we obtain for : 



-D(r) cos 



7rr- 



cos 



27rr- 



B 



£>(r) cos 



7rr- 



sm 



27rr- 



+ T COS 

4 



+ - sin 

4 



'"'■-B- + 4 



(50) 



(51) 



Here, -Fmax,-Pmin correspond to the maximum and minimu m cross-sectional areas of the FS, respectively. 
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